Scattering universality classes of side jump in anomalous Hall effect 
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The anomalous Hall conductivity has an important extrinsic contribution known as side jump 
contribution, which is independent of both scattering strength and disorder density. Nevertheless, we 
discover that side jump has strong dependence on the spin structure of the scattering potential. We 
propose three universality classes of scattering for the side jump contribution, having the characters 
of being spin-independent, spin-conserving and spin- flip respectively. For each individual class, the 
side jump contribution takes a different unique value. When two or more classes of scattering are 
present, the value of side jump is no longer fixed but varies as a function of their relative disorder 
strength. As system control parameter such as temperature changes, due to the competition between 
different classes of disorder scattering, the side jump Hall conductivity could flow from one class 
dominated limit to another class dominated limit. Our result indicates that magnon scattering plays 
a role distinct from normal impurity scattering and phonon scattering in the anomalous Hall effect 
because they belong to different scattering classes. 

PACS numbers: 72.10.-d,73.50.Bk,05.30.Fk,72.25.-b 



I. INTRODUCTION 

The anomalous Hall effect (AHE), in which a trans- 
verse voltage is induced by a longitudinal current flow 
in ferromagnetic materials, is one of the most intriguing 
effects in condensed matter physics. While it has been 
widely used experimentally as a standard technique for 
the characterization of ferromagnet materials, the theo- 
retical formulation of AHE proves to be complicated and 
is a subject full of controversial issues and conflicting re- 
sults.— In recent years, an important connection has been 
established between AHE and the Berry phase of Bloch 
electrons^— This triggers revived interest in this subject 
and is followed by extensive researches both theoretically 
and experimentally^— 

It is now generally accepted that spin-orbit coupling 
and spin splitting are two essential ingredients for AHE, 
and apart from an intrinsic contribution which is scatter- 
ing independent, there are also important extrinsic con- 
tributions to AHE due to disorder scattering. Based on 
their parametric dependence on the disorder density ndi SJ 
the extrinsic contributions can be collected into two sub- 
groups: the side jump contribution of order n[J is and the 
skew scattering contribution of order n~^~ s . 

The side jump contribution is of special interest in that 
it arises from scattering, but surprisingly it does not de- 
pend on either the scattering strength or the disorder 
density (we shall use the term "disorder strength" to 
stand for both scattering strength and disorder density). 
Theoretical calculations of simple model systems show 
that the side jump contribution is usually at least as im- 
portant as intrinsic contribution! 16 ' 17 i 22 ' 24 However, the 
good agreement between the intrinsic contribution cal- 
culated from first principles and the experimental results 
seems to indicate that the side jump contribution only 



plays a subdominant rolei 13 i 19 This remains as a puzzle 
that need to be resolved and is partly our motivation for 
the present work. 

Historically, the concept of side jump was first devised 
by Berger^ which refers to the coordinate shift of a 
wave-packet during an impurity scattering and this pro- 
cess leads to a contribution of order n^ is to the anomalous 
Hall conductivity. Recently, it has been found that be- 
sides this coordinate shift process, several other scatter- 
ing processes also generate contributions of order n® Us . M 
It should be noted that the term "side jump" used in 
the present paper includes all the scattering induced con- 



tributions of order 



Mis' 



not only the contribution from 



Berger's original side jump. 

In the study of physical systems, properties which are 
insensitive to detailed parameter values and system con- 
figurations but are only determined by the symmetry 
are especially interesting and important. It is helpful 
to define universality classes based on the behavior of 
these universal properties under certain imposed symme- 
try and study the generic properties of each class. Side 
jump can be regarded as a universal property for a disor- 
dered system in the sense that its value does not depend 
on the detailed disorder profile, but we shall see that it 
has sensitive dependence on the symmetry property of 
the scattering. Consequently, it is natural to define uni- 
versality classes of disorder scattering according to their 
side jump contributions and study the anomalous Hall 
response for each class. 

In this work, we propose three universality classes of 
disorder scattering, each has different structures in spin 
space. We find that: (1) for each individual class, the 
side jump contribution takes a distinct value indepen- 
dent of the detailed disorder profile. In particular, we 
show that magnon scattering plays a distinct role from 
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both impurity scattering and phonon scattering in AHE; 
(2) when several classes of scattering are present, side 
jump depends on their relative disorder strength and a 
sign change is possible as a result of their competition. 
Since in real physical system scattering processes of all 
the three classes exist, our finding indicates that a careful 
classification and analysis of different scattering processes 
is indispensible for an accurate account of AHE. 

This paper is organized as the following. First, in Sec. 
II, we propose and discuss the three scattering univer- 
sality classes of side jump. In Sec. Ill we demonstrate 
our ideas by a concrete analytical calculation of anoma- 
lous Hall conductivity of massive Dirac model for each 
universality class. In Sec. IV we discuss the important 
consequences of our result, especially about contribution 
from magnon scattering in AHE, and draw some final 
conclusions. 



II. UNIVERSALITY CLASSES OF DISORDER 
SCATTERING 

The general form of a random disorder potential for 
carriers with spin (or pseudospin) degrees of freedom can 
be written as 

Vdis(r) = [ y o( r - R >) + V ( r - ^) ■ °\ > 

i 

where Ri (i = 1,2,- ••) are the positions of randomly 
distributed scattering centers, a = (<J x ,a y ,a z ) is a vector 
with components of Pauli matrices and the hat means the 
quantity is a 2 x 2 matrix in spin space. We assume that 
the statistical average of the disorder potential is zero 
(any nonzero value only shifts the origin of total energy) 
and the second order spatial correlation only depends on 
the difference in positions, 

(V dis (r)) c = 0, (V dis (r)V dis (r')) c = B(r-r'), (2) 

where the angular bracket (• • ■ ) c denotes disorder aver- 
age. In order to discuss skew scattering which originates 
from higher order scattering processes, we allow a non- 
vanishing third order disorder correlation instead of re- 
quiring the disorder to be purely Gaussian. 

Time reversal symmetry has to be broken for the ap- 
pearance of AHE. 1 In ferromagnet, this is realized by 
the spontaneous magnetic ordering. We will be most in- 
terested in the configuration that the magnetization is 
perpendicular to the two-dimensional (2D) plane where 
the transport occurs, as is pertinent for most experimen- 
tal investigations. It is reasonable to assume that over 
disorder average the system is isotropic in the 2D plane 
with no preferred in-plane directions. With this sym- 
metry constraint, the total angular momentum (in the 
direction normal to the plane which we will refer to as 
z-axis) of the carriers are conserved on average. Due 
to spin-orbit coupling, the carrier's orbital motion which 
is tied to orbital angular momentum depends sensitively 



on the change of its spin angular momentum during a 
scattering. Based on this consideration, we propose the 
following three classes of disorder scattering which as we 
will see lead to different values of side jump contribution: 

Class A V = V°l, 

Class B V = V°a z , (3) 

Class C V = V°a±/V2, 

where V° denotes the orbital part of the scattering po- 
tential, and a± = a x ± ia y . Each class has a different 
action on the carrier's spin. Class A is isotropic in spin 
space. Class B, like Class A, conserve the z-component 
of the carrier spin but spin up and spin down carriers ex- 
perience different scattering potentials. Class C, unlike 
the first two classes, induces spin flips. The three classes 
as we discuss later represent a quite general classifica- 
tion scheme for real physical systems. This classification 
scheme is also evident if we consider the disorder correla- 
tion function under the in-plane rotational symmetry^ 

{v d 1 a (r)Vg(r')) c = (V„V )c% + (V X VM&M&*)# 

cx.—x,y 

(4) 

where i,j are spin indices. Since (V X V X ) C — (V y V y ) c due 
to the in-plane rotational symmetry, the last term is pro- 
portional to 

[{v+)ij{cr-)j l + {a-)ij{o + ) }i ] , (5) 

with each term being invariant under spin rotations 
around z-axis. The three terms in Eq.Q just correspond 
to the three classes we defined. 

Before proceeding, we point out an important differ- 
ence between Class C and Class A, B on their third order 
correlation functions. The Class C disorder can be ex- 
pressed as V(r)-cr where V is a random in-plane vector. 
Under in-plane rotational symmetry, V has no preferred 
direction therefore its third order correlation like (VVV) C 
must vanish. However for Class A or Class B, the third 
order correlation is not dictated by this symmetry con- 
straint hence does not necessarily vanish. This difference 
will be reflected in the skew scattering contribution to 
the AHE. 

The transverse motion of carriers in AHE is a result 
of spin-orbit coupling. In our classification scheme, each 
class of scattering has different effects on the carrier spin, 
hence will also have different effects on the carrier orbits. 
This is the underlying reason for their distinct contri- 
butions to the AHE and especially the side jump part. 
In the following section, we demonstrate this idea by a 
concrete model calculation. 
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III. AHE OF MASSIVE DIRAC MODEL 
A. Model and Approach 

To demonstrate the rationale of our classification 
scheme, we calculate the anomalous Hall conductivity for 
the massive Dirac model. This model is usually consid- 
ered as the minimal model for AHE. 1 The model Hamil- 
tonian reads (we set H = 1 and assume A > in the 
following calculations) 



H = v(k x a x 



ky&y 



Aa z 



(6) 



where spin-orbit coupling is contained in the first term 
with v being the coupling constant, and the last term 
breaks the time reversal symmetry and is also responsi- 
ble for the finite electron mass at the band edge. This 
model captures interesting physics near a generic band 
anti-crossing point due to spin-orbit coupling. 
The eigenstates of the system are 



1 



ik-r I 



with the corresponding energy eigenvalues 



e ± (k) — ±y/ (vk) 2 + A 2 , 



(7) 



(8) 



where ± labels the upper and lower band respectively, 
A is the system size, and \u k ) is the spin part of the 
eigenstate which can be written as 



cos 2 
= 1 sinfe 1 * 1 



sin ■ 



, (9) 



where 8 and <f> are the spherical angles of the vector 
(vk x ,vk y , A) such that 



A 



y/(vk) 2 + A 2 ' 
tan <j) = k y /k x . 



sin 9 



vk 



^(vk) 2 + A 2 (10) 



Due to spin-orbit coupling, the spin state is a function 
of the momentum k. The energy spectrum consists of 
two anti-crossing bands with a band gap of 2A. From 
the dispersion relation Eq. © , the geometry of the bands 
can be termed as a "Dirac hyperboloid" (of two sheets). 

To calculate the anomalous Hall conductivity, we 
follow Sinitsyn et al£^ by using the Kubo-Streda 
formalismj&^ In this approach, the Hall conductivity 
can be separated into two parts in the weak scattering 
regime, a xy = a xy + a x l y , where a xy is a Fermi surface 
contribution which includes all the important scattering 
contributions, and a l xy is a Fermi sea contribution for 
which wc only need to retain the scattering-free compo- 
nent^ In the following we consider that the system is 
electron doped with Fermi energy sp > A and due to 
particle-hole symmetry, the results can be easily general- 
ized to the hole-doped case. We assume that the system 



is in the weak scattering regime, i.e. kpl ^> 1 where kp 
is the Fermi wave vector and I is the electron mean free 
path. It has been found that a l xy vanishes^ and the task 
gets reduced to the evaluation of a xy which is given by 
the following expression 



2nA 



TY(v x G R (e F )v y G A (e F ) 



(11) 



where G R and G A are the retarded and advanced Green's 
functions respectively, v x and v y are the velocity opera- 
tors, and the trace is taken over both momentum and 
spin spaces. In weak scattering regime, the calculation is 
performed perturbatively in the small parameter 1/ (kpl). 

In our model, we consider the scattering processes to 
be quasi- elastic, hence the disorder lines in Feynman dia- 
grams carry no energy arguments. This serves as a good 
approximation for the scattering by collective excitations 
such as phonons or magnons as long as energy of collec- 
tive excitation involved in the scattering is much less than 
the Fermi energy. Since the typical energy scale of excita- 
tions is fegT, this condition is satisfied for temperatures 
with ksT <C ep. Furthermore, for massless excitation 
with a spectrum u)(k) — v q k (v q is a constant sound 
speed), quasi-elastic approximation is justified even at 
higher temperatures if the quasi-particle speed v q is much 
less than Vp, viz the band velocity at Fermi level. For 
massive quasi-particle excitations, its validity can be jus- 
tified if the quasi-particle mass is much larger than the 
electron effective mass. 

In the following, we calculate the Hall conductivity for 
each individual class, or equivalent ly when one class of 
scattering is dominant. The evaluation of the conductiv- 
ity follows standard procedures, and the relevant Feyn- 
man diagrams under self-consistent non-cross approxima- 
tion have been identified before^ so we do not elaborate 
here. The results are listed below. 



B. Intrinsic Contribution 

The intrinsic contribution of AHE is a property purely 
of the spin-orbit coupled band structure. It was first 
proposed by Karplus and Luttinger^ and recently its 
connection with the Berry phase of Bloch electrons is 
established^^ It is now understood that spin-orbit cou- 
pled bands usually possess effective magnetic fields in 
momentum space known as Berry curvatures^ which de- 
flect carriers in the transverse directions. The intrinsic 
contribution of anomalous Hall conductivity equals the 
integration of Berry curvatures of all the occupied states. 
Because it does not depend on scattering, intrinsic contri- 
bution is the same for all the three universality classes. 
In Kubo-Streda formalism, intrinsic contribution is the 
sum of the scattering- free part of a 1 and a x l y . 

For electron doped case, we can separate the intrinsic 



Hall conductivity a X y into two parts, 



int(v) 



int(c) 



(12) 



4 



T 



where Of is the spherical angle 9 at the Fermi surface 
when k — hp. 



(a) (b) 




FIG. 1: (a)-(d) are the four conductivity diagrams of side 
jump that correspond to the distribution function correction. 
+ and — represent the upper and lower band respectively. T 
stands for the renormalized velocity vertex which is dressed 
by a ladder like diagram as shown in (e). 

where Oxy^ is the contribution from all the completely 
occupied valence bands below the Fermi surface and 
(jxy^ 1 is the contribution from the partially filled con- 
duction band where the Fermi surface lies in. The con- 
tribution from completely filled bands Cxy^ must be a 
topologically quantized value Ce 2 /(27r) with C being an 
integer called the first Chern number^ The lower band 
of the massive Dirac model has a contribution of — e 2 /47r. 
This is not a contradiction because the Dirac band is not 
bounded. For any real physical system, the evaluation 
of C must go beyond the low energy effective model and 
require complete information of the entire Fermi sea. On 
the contrary, the contribution a^y^ from the partially 
filled conduction band can be regarded as a Fermi surface 
property^ and is captured within the effective model, 

< (C) = (13) 



C. Side Jump Contribution 

Now let's focus on the side jump contribution which 
is the central quantity we are interested in. For each 
individual class, it is independent of disorder density ridis 
and scattering strength V° . It can be expressed in terms 
of 6p and a set of scattering times defined on the Fermi 
surface. For notational convenience, we define 

L = 2nn dis J ^\ V £, k \ 2 cos\4~^)S(e F -e+), (14) 

where V£, k — (k'\V°\k) is the matrix element of the or- 
bital part of the scattering potential in momentum space 
and i = 0, 1, 2, • • • is an integer. 

In the semiclassical picture, the side jump we defined 
here consists of three components: a contribution from 
the coordinate shift (the original Berger's side jump), a 
contribution from a correction of distribution function 
(called the anomalous distribution), and a contribution 
from some higher order scattering processes (called the 
intrinsic skew scattering). The first two components are 
shown to be equali 2 ^ In the Kubo-Streda approach, Fig.Q] 
shows a set of diagrams that contributes to the side jump 
in the chiral (eigenstate) basis. These correspond to the 
contribution from the anomalous distribution function 
correction, i.e. the second component above. The dia- 
grams corresponding to the contribution from coordinate 
shift can be obtained by simply exchanging the subscripts 
x and y in Fig. [T] and further making a 180° rotation (i.e. 
exchanging G R and G A ). The resulting contribution to 
Hall conductivity from these two components for each 
scattering class is 



Class A 
Class B 
Class C 

We observe that different scattering class contributes very differently to the Hall conductivity. In the diagramatic 
approach, this difference originates from the different k dependence at scattering vertices, which in turn results from 
their different spin structures. It should be noted that the vanishing value of class B is not a general feature but 
rather depends on the specific model we considered here4£ 



sin 



T s K a ) 

xy 2 7 t(1 + cos 2 ^)t - 1 

T s K a ) _ n 

'xy 

„2 

7 s j(a) 
'xy 



• cos 8p (tq -t x ) 



2 cos 2 9ft x 



sin 2 9ft 2 1 



(15) 



47T 



cos 9 1 



5 




(e) (f) 

FIG. 2: The diagrams corresponding to the part of side jump from fourth order scattering process (intrinsic skew scattering). 



The third component (intrinsic skew scattering) results from certain fourth order scattering processes. The corre- 
sponding diagrams are shown in Fig. [2] and the results are 



Class A : 

Class B : 
Class C : 



T sj(b) 



T sj(b) 
'xy 



e 2 sin 4 9p cos 9f(t 



2t; 



4lT [(1 + cos 2 Op)^ 1 - 2 cos 2 9ft- 1 - sin 2 f t 2 -1 ] ' 



sin 4 9p cos #f(t ( 



-l 
o 



.-1X2 



(16) 



47r [(1 + cos 2 0f)t o -1 - 2t^ + sin 2 e F T^] ' 



CT sj(b) = 

w xy w ' 



The a%y contribution vanishes for class C is a very gen- 
eral result because each such diagram contains a factor 
of the form J sin <j>d(f> from the momentum integral at the 
velocity vertex which suppresses the intrinsic skew scat- 
tering process. 

The total side jump contribution to anomalous Hall 
conductivity is given by a^ y — axy + Oxy . It is 
clear that each class has a distinct side jump contribu- 
tion. Scattering rates with the same power appear in 
both nominator and denominator of the expressions in 
Eqs. p5ll6[) , hence the results are independent of disor- 
der density and scattering strength. The dependence on 
the band parameters such as A and Ef are also different 
for different classes. Furthermore, it should be noted that 
different classes can have side jump contribution with dif- 
ferent signs, as shown here between Class A and the other 
two classes. 



tribution for each scattering class. In the semiclassi- 
cal picture, skew scattering contribution comes from the 
asymmetric part of the scattering rates for higher order 
scattering processes. The leading contribution is related 
to the third order disorder correlation and has a nZ 1 




(a) (b) 



FIG. 3: The two diagrams corresponding to skew scattering 
contribution. 



D. Skew Scattering Contribution 



Although our focus is the side jump contribution, to 
be complete, we also calculated the skew scattering con- 



dependence 4i The corresponding Fcynman diagrams are 
shown in Fig. [3] and the skew scattering contribution to 
Hall conductivity for each class is given by 
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Class A : 

Class B : 
Class C : 



where 



— = 2ne F n d [ s 



sk 
'xy 



sk 



sin Of cos Oft, 



.sk 



27r [(1 + cos 2 Op)^ 1 - 2 cos 2 F Tf 1 - sin 2 9 f t^ x ) 2 ' 



sin 4 9 f tJ 



(17) 



2tt 



[(1 



cos- 1 w F it, 



- 2t _1 
o ZT i 



= 0. 



cPk' 



d 2 k 

W) 

sm{<t> — 



sm 



sin 2 6pT 2 1 ] 2 



5<V£„V£ k ,,V£,J c siii(0'-0) 



- 0') + sm(^' 
I 



(18) 



S(e F - e%,)6(s F - e£„). 



Note that the factor l/r 2 k is proportional to ndis, hence Class C: 
(T* is of order nj is . Because the third order disorder 
correlation vanishes for Class C as we discussed in Sec. II, 
V T sk c = ^ an d the skew scattering process is forbidden 
for this type of disorder scattering. Therefore, we see an 
important qualitative difference of the anomalous Hall 
response for Class C and the other two classes: in the 
weak scattering regime, the leading contribution for Class 
A and Class B is the skew scattering of order n~f , but 
for Class C the leading contribution is the intrinsic plus 
side jump which are of order n^ is . 



E. Total Hall Conductivity of Order n° dis 

The above results are valid for disorder scattering with 
general orbital part. If we consider the simple white noise 
(short range) disorders, the results for Hall conductivity 
are greatly simplified. In this case, we have for each class 



1 o 11 

Tl T 2 2t ' 



(19) 



Then the total Hall conductivity can be written as (o^ 
is not included here, as discussed in Sec. III. B): 
Class A: 



47T 



T — cos 6 F 



e 2 sin 2 F cos Op 
it 1 + 3 cos 2 Op 



3e 2 sin Op cos Op 



2e 2 



sm VpcosOp t. 



sk 



4tt (1 + 3cos 2 F ) 2 7T (1 + 3cos 2 (9f) 2 T - 



(20) 



Class B: 



Oxy =-j— (1 - C0S6> F ) + 



47T 



e 2 sin 4 Op cos F 
4^ (3 + cos 2 Op) 2 



+ 



2e 2 



sin 4 Of 



(21) 



'sk 



7T (3- 



?f) 2 



e 

47r' 



(22) 



The expression in Eg. ((20]) recovers the result obtain by 
Sinitsyn et al— Clearly, each universality class has its 
distinct extrinsic Hall conductivity and different func- 
tional dependence on the system parameters such as Of. 
Note that for class C, the side jump contribution can- 
cels with the part of intrinsic contribution which de- 
pends on the Fermi energy, making the final result a 
constant and the value is the same as the intrinsic con- 
tribution for a completely filled conduction band (i.e. in 
the limit Of —> n/2). We have checked that this interest- 
ing cancelation occurs for a generic class of Hamiltonians 
rl = vk n [cos(n(j)k)(j x + wa.(n4>k)cFy\ + &-<J z where n is an 
integer. 

Here we are most interested in the part of Hall conduc- 
tivity a® y that is of order n dis . This includes the intrinsic 
contribution and the side jump contribution, i.e. 



_int 



(23) 



In Fig. |4j for each scattering class (with white noise spa- 
tial correlation) , we plot <r° y as a function of Fermi en- 
ergy £f for the massive Dirac model. Observe that cr® y 
for class C takes a constant value e 2 /47r which is inde- 
pendent of Fermi energy and the curves for both class 
A and class B approach this constant value asymptoti- 
cally as ef — > oo. For class A, the extrinsic contribu- 
tion has opposite sign as compared with the the intrinsic 
contribution. Because of this sign difference, the Hall 
conductivity for class A takes negative values for Fermi 
energies below £p s=s 7.3A. This behavior differs from 
that for class B and class C whose extrinsic contribu- 
tions have the same sign as the intrinsic contribution, so 
their overall <r° y 's are positive. This shows that the Hall 
conductivity sensitively depends on the class of disorder 
scattering. 
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FIG. 4: (color online). a° plotted as a function of the Fermi 
energy ef for each of the three universality classes. a° y is 
measured in units of e 2 /(47r), and ef is measured in units of 
A which is half of the band gap. 

F. Competition between Classes 

In the presence of two or more classes of scattering, 
there will be a competition between them in the anoma- 



lous Hall response. The resulting side jump contribution 
takes the following generic form: 



.sj E a a a T * 1 C a pT a 1 T fs 

xy T—y i — 1 ' i — 1 — 

2^ a °a r a a a0 T a Tg 



where the is the scattering time defined for each 

class of scattering involved, a a , b a , c a /3 and d Q( g are the 
(disorder independent) coefficients which depend only on 
system intrinsic parameters such as 8p in the present 
model. 

As an example, let's consider the competition between 
Class A and Class C. The calculation is tedious but 
straightforward. The resulting total Hall conductivity 
can be expressed as 



e 2 , A s e 2 

T M , = (1 — COS Op) 

y 4tt V b > 7r(l + 3cos 2 F 



sin 2 Of cos # f (1 — C) 



sin 4 F cos6> F (§ - C + 2r?) 



4sin^6» F C * [(l + 3cos 2 6» F )+4sin 2 6» F C]' 



(25) 



where the parameter C, defined as £ = (r^i — t^ c )/ (t^ — 
t^2) is a measure of the relative disorder strength of 
the two classes, and r\ = t~^ a /{tq^ — t^ 1 ) 2 is a factor 
for skew scattering contribution from Class A, here 
stands for the scattering time Tj defined in Eq. (fT4"| for 
Class A scattering and nc is similarly defined. The first 
term above is the intrinsic contribution, and the remain- 
ing two terms (with 77 = 0) are the side jump contri- 
bution. Observe that in the limit ( — > or ( — > +00, 
Eg. (1251) recovers our previous results in Eg. (l20[) and 
Eq. (|2"21 , and the value of Hall conductivity varies contin- 
uously as £ changes between these two limits. This shows 
that the value of side jump is no longer independent of 
disorder strength but can vary as a result of competition 
between different scattering classes. 

In Fig. [5l we plot the Hall conductivity cr° y (by set- 
ting 77 = 0) as a function of the Fermi energy for differ- 
ent values of £. As C increases from zero, the curve of 
cr° y is shifted upward from the Class A dominated case 
due to the increasing contribution from Class C scatter- 
ing, and finally approaching the value e 2 /(47r) for the 
Class C dominated case. This competition behavior is 
more clearly seen in Fig. [6j where a® y is plotted at three 
fixed Fermi energies as a function of C- We see that as 
£ increases, (T® y increases monotonically. In the energy 
range e < 7.3A, there is a sign change of a® y during this 
crossover. 
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FIG. 5: (color online). a xy plotted as a function of the Fermi 
energy ef for fixed values of C,. a^ y is measured in units of 
e 2 /(47r), and ef is measured in units of A. 



IV. DISCUSSION AND CONCLUSION 

As demonstrated in Sec. Ill, different scattering class 
has its own distinct contributions to AHE. This suggests 
that for the study of AHE in real materials, competing 
scattering processes belonging to different classes need to 
be handled carefully. 

In ferromagnetic materials, normal (non-magnetic) im- 
purity scattering and phonon scattering belong to Class 
A since they are isotropic in spin space. Most of the 
previous studies on extrinsic AHE are focused on this 
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FIG. 6: (color online). a® y versus the ratio of scattering rates 
£ for fixed values of Fermi energy ef- &xy is measured in 
units of e 2 /(47r). The plot shows the crossover from Class A 
dominated regime to Class C dominated regime as £ increases. 

class of scattering and indeed it has been found that 
the electron-phonon scattering has similar contribution 
as normal impurity scattering (although there is no skew 
scattering due to conservation of phonon population in 
steady state) which is consistent with our theory. Mag- 
netic impurities with spin directions oriented along the 
average magnetization is of Class B and they should gen- 
erate a contribution different from that of normal impu- 
rities. This has also been observed in the study of dilute 
magnetic semiconductors ! 31 ' 49 

Scattering processes of Class C also exist. For exam- 
ple, magnetic impurities with random in-plane magnetic 
orientation are of this class. Moreover the scattering of 
electron by magnons also belongs to Class C. To see this 
more clearly, let us consider the coupling between con- 
duction electron spin <x and the local spin S (within an 
s-d model approach), 

Hint = -J fdr [<r(r) • S(r)} 

j r (26) 

= -- dr (a+S- + a_S+ + 2& Z S Z ) , 

where J is the exchange coupling constant. The last term 
2a z S z describes the Zeeman splitting (which has already 
been included in non-perturbed part of the Hamiltonian) , 
whereas the first two terms describe the scattering by 
magnetic excitations. There is a transfer of spin angular 
momentum between conduction electrons and local spins 
during this process, hence such scattering is of Class C 
in our classification. 

For real material samples that are studied experimen- 
tally, all the three classes of scattering are present at 
finite temperature. At low temperature, the scattering 



by normal impurity usually dominates. With increasing 
temperature, electron-magnon scattering becomes more 
important and can compete with normal impurity scat- 
tering and phonon scattering. Especially for materials 
with a high Debye temperature and a low Curie temper- 
ature, we can conceive a situation in which Class A and 
Class C scattering compete as depicted in Sec. III.F. 
Then the value of the side jump Hall conductivity would 
flow between the two limiting values as a function of tem- 
perature. 

Finally we point out that the concept of "spin" in 
our discussion can be very general, corresponding to any 
discrete degrees of freedom (sometimes called a "pseu- 
dospin"). For example, in a bipartite lattice (such 
as graphene), the sublattice degree of freedom can be 
treated as pseudospin. Anomalous valley Hall transport 
occurs in graphene when there is sublattice symmetry 
breaking in the system^ For bilayer systems, it is the 
layer index that plays the role of pseudospin. Then scat- 
tering processes can be classified according to their effects 
on the pseudospin. For example, inter-sublattice scat- 
tering in a bipartite lattice or interlayer scattering in a 
bilayer system would both belong to Class C. In general, 
our results indicate that a careful analysis of various scat- 
tering processes according to their pseudospin structures 
is indispensable for the study of AHE in these systems. 

In summary, we have shown that the extrinsic part 
of the anomalous Hall conductivity has a strong depen- 
dence on the spin structure of the disorder scattering. 
We propose three universality classes of scattering ac- 
cording to their side jump contribution to the anomalous 
Hall conductivity. Each class has its distinct value of 
side jump. When two or more classes of scattering are 
competing, the side jump contribution is determined 
by their relative disorder strength. Various scattering 
processes in real physical systems can be classified into 
these three classes. In particular, we demonstrate that 
magnon scattering has distinct side jump contribution 
from normal impurity scattering and phonon scattering 
and the value of side jump contribution could change as 
system control parameter (such as temperature) varies. 
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